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Abstract 

In this paper we determine the number of the meaningful compositions of higher 
order of the differential operations and Gateaux directional derivative. 

1 The compositions of the differential operations 
of the space R 3 

In the real three-dimensional space IR 3 we consider the following sets: 

A = {/:R 3 — ^M|/gC ,00 (M 3 )} and A ± = {f.R 3 — >R 3 | fe G°°(M 3 )}. (1) 

Then, over the sets Ao and Ai in the vector analysis, there are m = 3 differential 
operations of the first-order: 

curl/ = V 2 J= %_|A\ ; Ai _ Ai (2) 

\ox2 0x3 0x3 ox\ ax\ 0x2 J 



dxi 8x2 dx 
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Let us present the number of the meaningful compositions of higher order over the 
set A3 = {Vi,V2,Vs}. As a well-known fact, there are m = 5 compositions of 
the second-order: 

A/ = divgrad/ = V30VJ, 
curl curl / = V2 V2 /, 

graddiv/= VxoVa/, (3) 
curl grad / = V2 Vi / = 0, 

divcurl/ = V 3 o V 2 /= 0. 
Malesevic [2] proved that there are m = 8 compositions of the third-order: 

grad div grad / = Vi o V3 o Vi /, 
curl curl curl / = V2 V2 V2 /, 
div grad div/ = V3 o Vi o V3/, 
curl curl grad / = V2 V2 Vi / = 0, 

(4) 

div curl grad / = V3 o V2 Vi / = 0, 
div curl curl / = V3 o V2 V2 / = 0, 
grad div curl / = Vi o V3 o V2 / = 0, 
curl grad div / = V2 Vi o V3 / = 0. 

If we denote by f(k) the number of compositions of the /c th -order, then Malesevic [3] 
proved: 

f (fe) = F k+3 , (5) 

where is k th Fibonacci number. 

2 The compositions of the differential operations 
and Gateaux directional derivative on the space R 3 

Let / e A be a scalar function and e = (ei, e^) 6 I 3 be a unit vector. Thus, the 
Gateaux directional derivative in direction e is defined by [1, p. 71]: 

df df df 
div s f = V / = Vxf ■ e = -Le ± + ^-e 2 + e 3 : A — . A . (6) 

OX 1 OX2 OX3 
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Let us determine the number of the meaningful compositions of higher order over 
set £> 3 = {Vo, Vi, V2, V3}. There exist m = 8 compositions of the second-order: 

dir e - dir- / = Vo o Vo / = Vi ( Vi/ • e ) • e, 

grad dir- / = Vi o V / = Vi ( VJ ■ e ) , 

A/ = divgrad/ = V 3 oVi/, 

curl curl / = V2 V2 /, 

dir g div/= V oV 3 /= (V10V3/) -e, 

grad div / = Vi o V 3 /, 

curl grad / = V2 Vi / = 0, 

div curl/ = V 3 oV 2 / = 0; 
that is, there exist m = 16 compositions of the third-order: 



Hit Hit Hit t — 

air- Qll^. Qllg. J - 


= Vo 


Vo 


0V0/, 




grad du% dir^ / = 


= Vj 


V c 


0V0/, 




div grad di% / = 


= v 3 


Vi 


0V0/, 




dir- div grad / = 


= Vo 


v 3 






grad div grad / = 


= Vi O V; 






curl curl curl / = 


= v 2 


v 2 


0V2/, 




dir- dir^. div / = 


V c 


> V < 


>V 3 /, 




grad dir ? div / = 


= Vi 


Vo 


0V3/, 




div grad div / = 


V30V10 v 3 /, 




curl grad dh% / = 


= Vj 


,0 V] 


.oVo/ = 


= 0, 


curl curl grad / : 


= V 


,0 v 


»oVJ : 


= 0, 


div curl grad / = 


= v 3 


v 2 


oVi/ = 


= 0, 


div curl curl / = 


v 3 « 


3 V 2 ' 


3V 2 / = 


0. 


dir- div curl / = 


V V 3 < 


>v 2 / = 


0. 


grad div curl / = 


= Vi 


v 3 


oV 3 /= 


= 0, 


curl grad div / = 


= v 2 


Vi 


oV 3 / = 


= 0. 



Using the method from the paper [3] let us define a binary relation o u to be in 
composition" : Vj cr Vj = T iff the composition Vj o V, is meaningful. Thus, Cayley 
table of the relation a is determined with 



a 


V 


Vi 


v 2 


v 3 


V 


T 


T 


JL 


JL 


Vi 


JL 


1 


T 


T 


v 2 


JL 


1 


T 


T 


v 3 


T 


T 


1 


J_ 



Let us form the graph according to the following rule: if V« a Vj = T let vertex Vj 
be under vertex Vi and let there exist an edge from the vertex Vi to the vertex Vj. 
Further on, let us denote by V_i nowhere-defined function i?, where domain and range 
are the empty sets [2]. We shall define V_i a Vj = T (i = 0, 1, 2, 3, 4). For the set 
£> 3 U {V_i} the graph of the walks, determined previously, is a tree with the root in 
the vertex V_i. 




g(o) 
g(i) 

g(2) 
g(3) 
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Fig. 1 

Let g(/c) be the number of the meaningful compositions of the /c th -order of the functions 
from B3. Let gj(fc) be the number of the meaningful compositions of the /c th -order 
beginning from the left by Vi. Then g(k) = g (k) + g 1 (fc) + g 2 (^) + g 3 (^)- Based on 
the partial self similarity of the tree (Fig. 1) we get equalities 



(10) 



&>(*) 


= &(*-!) 


+"gi(*-l 


g#) 


= &,(*-!) 


^g 3 (fc-l 


&(*) 


= &(*-!) 


f g 3 (A;-l 


gs(*) 


= &(*"!) 


H-gi(*-l 



Hence, a recurrence for g(/c) can be derived as follows: 

g(A;) = 2g(*-l). 
Based on the initial value g(l) = 4, we can conclude: 



g(*) 



(12) 
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3 The compositions of the differential operations 
of the space R n 

Let us present the number of the meaningful compositions of differential operations 
in the vector analysis of the space MJ 1 , where differential operations V r (r = 1, . . . , n) 
are defined over non-empty corresponding sets A s (s = 1, . . . , m and m— \n/2\ , n > 3) 
according to the papers [3], [4]: 

A n (n = 2m + l): Vi : A ^A 1 



An (n = 2m): Vi : A ^Ai 
V 2 : Ai^A 2 

Vi : A*_i->Ai 

Vm • A m _i > A m 
V m +i : A m — >A m -i 

V n -j : Aj +1 — ► Aj 



V 2 : Ax^A 2 

Vi : Ai_x^Ai 

V m . A m „i > A m 

V m+ i : A m — >A m (13) 

V m + 2 : A m >A m _i 

V n _j : Aj + i — > Aj 



V n _i : A 2 ^Ai 
V n : A X ^A , 



V n _i : A 2 ^Ai 
V n : Ax^A . 

Let us define higher order differential operations as the meaningful compositions of 
higher order of differential operations from the set A n = {Vi, . . . , V„}. The number 
of the higher order differential operations is given according to the paper [3]. Let us 
define a binary relation p 11 to be in composition": VipVj = T iff the composition 
Vj o Vi is meaningful. Thus, Cayley table of the relation p is determined with 

V iP V j =\ . (14) 

_L , otherwise. 

Let us form the adjacency matrix A = [ay] G { 0, l} nxn associated with the graph, 
which is determined by the relation p. Thus, according to the paper [4], the following 
statement is true. 

Theorem 3.1. Let P n (X) = |A — AI| = aoX n + aiA* 1-1 + ■■■ + «„, be the characteristic 
polynomial of the matrix A and v n — [ 1 . . . 1 ]i X n- we denote by f(k) the number of 
the k th -order differential operations, then the following formulas are true: 

f{k) = v n -k k ~ 1 -v T n (15) 

and 

a f (k) + aif (k - 1) H h a n f (k - n) = (k>n). (16) 



5 



Lemma 3.2. Let P n (X) be the characteristic polynomial of the matrix A. Then the 
following recurrence is true: 



P n (A) = A 2 (P n _ 2 (A)-P n _ 4 (A)). 



(17) 



Lemma 3.3. Let P n (A) be the characteristic polynomial of the matrix A. Then it has 
the following explicit representation: 



,fe_l/ --k+2 



k=l 



-ir 



fc-i 



A 



— 2A+2 



lH 1 i+ 2 ( , n+3 v / 



n+3 
~~ 2 
fc-2 



A A 



n-2fc+2 



n = 2m; 



n = 2m + l. 



(18) 



The number of the higher order differential operations is determined by corresponding 
recurrence, which for dimension n = 3,4,5,...,10,we refer according to [3]: 



Dimension: 


Recurrence for the number of the /c th -order differential operations: 


n = 3 


f(fc) = f(fc-l) + f(fc-2) 


n = 4 


f (fc) = 2f (fc - 2) 


n = 5 


f (fc) = f (fc - 1) + 2f (fc - 2) - f (fc - 3) 


n = 6 


f(jfe) =3f(fc-2)-f(fc-4) 


n= 7 


f (fc) = f (fc - 1) + 3f (fc - 2) - 2f (fc - 3) - f (fc - 4) 


71= 8 


f(fc) = 4f (fc - 2) -3f(fc-4) 


71 = 9 


f (fc) = f (fc - 1) + 4f (fc - 2) - 3f (fc - 3) - 3f (fc - 4) + f (fc - 5) 


n = 10 


f (fc) = 5f (fc - 2) - 6f (fc - 4) + f (fc - 6) 



For considered dimensions n — 3, 4, 5, . . . , 10, the values of the function f (fc), for small 
values of the argument k, are given in the database of integer sequences [6] as sequences 
A020701 {n = 3), A090989 (n = 4), A090990 (n = 5), A090991 (n = 6), A090992 
(n = 7), A090993 (n = 8), A090994 (n = 9), A090995 (n = 10), respectively. 



4 The compositions of the differential operations 
and Gateaux directional derivative of the space 

Let / £ be a scalar function and e— (ei, . . . , e n ) G M n be a unit vector. Thus, the 
Gateaux directional derivative in direction e is defined by [1, p. 71]: 



8f 



fc=l 



An 



(19) 
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Let us extend the set of differential operations A n = {Vi,...,V n } with Gateaux 
directional derivational to the set B n = A n U {Vq} = {Vq, Vi, . . . , V n }: 



B n (n=2m): V : A ^A 
Vi : A ^A! 
V 2 : A X ^A 2 

Vi : Ai_i^Ai 



B n (n=2m+l): V : A ^A 
Vi : Ao^Ai 
V 2 : A!^A 2 

Vi : Ai„!^Ai 

Vm • A m _l > A r 



V m : A m _i — > A m (20) 



V m +i . A m ► A m _i 

V n _j : Aj + i^Aj 

V n -i : A 2 ^A X 
V n : A^Ao, 



V m + 2 : A m — ^A m _i 

V n _j : Aj +1 ^Aj 

V n _i : A 2 ^A X 
V n : Ai^A . 



Let us define higher order differential operations with Gateaux derivative as the mean- 
ingful compositions of higher order of the functions from the set B n = { Vo, Vi, . . . , V n }. 
We determine the number of the higher order differential operations with Gateaux 
derivative by defining a binary relation a "to be in composition" : 

( T , (z = 0Aj = 0) V (i = nAj=0) V (j = i + l) V (z+j = n + l); 
J_ , otherwise. 

Let us form the adjacency matrix B = [bij] E { 0, i}( n+1 ) xn associated with the graph, 
which is determined by relation a. Thus, analogously to the paper [4], the following 
statement is true. 

Theorem 4.1. Let Q n (X) = |B-AI| = f3 X n+1 + f3 1 X n H h/3 n +i be the characteristic 

polynomial of the matrix B andv n+ i = [1 . . . l]i x (n+i)- If we denote by g(k) the num- 
ber of the k th -order differential operations with Gateaux derivative, then the following 
formulas are true: 

g(k) = v n+l • B*" 1 ■ vl +l (22) 

and 

(3 g(k) + fcg(k - 1) + ■ • ■ + (3 n +ig(k - (n + 1)) = (k>n+l). (23) 

Lemma 4.2. Let Q n (X) and P n {X) be the characteristic polynomials of the matrices 
B and A respectively. Then the following equality is true: 

Q n (X) = X 2 P n - 2 (X)-XP n (X). (24) 
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Proof. Let us determine the characteristic polynomial Q n (X) = |B — AI| by 



Qn(X) 



A 


1 





. 


. 














-A 


1 


. 


. 








1 








-A 


1 . 


. 





1 














1 . 


. 


-A 


1 











1 


. 


. 





-A 


1 


1 


1 





. 


. 








-A 



Expanding the determinant Q n (X) by the first column we have 

Q n (X) = (1 - A)P n (A) + (-l) n+2 D n (A), 

where is 



Ai(A) 



1 








.. 














A 


1 





.. 











1 





-A 


1 


.. 








1 














1 . . 


. -A 


1 














1 


.. 





-A 


1 








1 





. . 








-A 


1 



(25) 



(26) 



(27) 



Let us expand the determinant D n (X) by the first row and then, in the next step, 
let us multiply the first row by —1 and add it to the last row. Then, we obtain the 
determinant of order n — 1 : 



1 








.. 











1 


A 


1 





.. 








1 








-A 


1 


.. 





1 














1 


.. 


. -A 


1 











1 





.. 





-A 


1 














.. 








-A 






D n (X) 



Expanding the previous determinant by the last column we have 



Dn(X) = (-I)" 



-A 





1 
-A 






1 

-A 1 



1 

1 




1 





-A 





1 
-A 1 

-A 



(28) 



(29) 
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If we expand the previous determinant by the last row, and if we expand the obtained 
determinant by the first column, we have the determinant of order n — 4 : 



D n (X) = (-1)"A 2 



-A 







1 

-A 1 



1 





-A 





In other words 



D n (X) = (-l)"A 2 P n _ 4 (A). 
^From equalities (31) and (26) there follows: 

Q n (A) = (l-A)P n (A) + A 2 iV 4 (A). 

On the basis of Lemma 3.2. the following equality is true: 

Q n (A) = A 2 P n _ 2 (A)-AP n (A). ■ 



1 

1 o 





1 
-A 1 

-A 



(30) 



(31) 
(32) 

(33) 



Lemma 4.3. Let Q n (X) be the characteristic polynomial of the matrix B. Then the 
following recurrence is true: 

Qn(A) = A 2 (Q n „ 2 (A)-Q n _ 4 (A)). (34) 

Proof. On the basis of Lemma 3.2. and Lemma 4.2. there follows the statement. ■ 



Lemma 4.4. Let Q n (X) be the characteristic polynomial of the matrix B. Then it has 
the following explicit representation: 



LiJ+i 



Qn{\) = { 



(X-2)J2 (— l) fe 1 ( I " J A" 
fc=i V 1 / 

f2 



L ni3 J+2 



k=l 



-2k+2 



n-2fc+3 



n = 2m + l; 



, n = 2m. 



(35) 



Proof. On the basis of Lemma 3.3 and Lemma 4.2. there follows the statement. 
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The number of the higher order differential operations with Gateaux derivative is 
determined by corresponding recurrences, which for dimension n — 3, 4, 5, . . . , 10 , we 
can get by the means of [5]: 



Dimension: 


Recurrence for the num. of the /c th -order cliff, operations with Gateaux derivative: 


n = 3 


g(k) = 2g(fc - 1) 


n = 4 


g(fc)=g(fc-l) + 2g(A;-2)-g(fc-3) 


n = 5 


g(fc) = 2g(fc - 1) + g(fc - 2) - 2g(fc - 3) 


n = 6 


g(fc) = g(fc - 1) + 3g(fc - 2) - 2g(fc - 3) - g(fc - 4) 


n = 7 


g(fc) = 2g(fc - 1) + 2g(fc - 2) - 4g(fc - 3) 


n = 8 


g(fc) = g(fc - 1) + 4g(fc - 2) - 3g(fc - 3) - 3g(fc - 4) + g(fc - 5) 


n= 9 


g(fc) = 2g(fc - 1) + 3g(fc - 2) - 6g(fc - 3) - g(fc - 4) + 2g(fc - 5) 


n =10 


g(fc) = g(fc - 1) + 5g(fc - 2) - 4g(fc - 3) - 6g(fc - 4) + 3g(fc - 5) + g(fc - 6) 



For considered dimensions n — 3, 4, 5, . . . , 10, the values of the function g(fc), for small 
values of the argument k, are given in the database of integer sequences [6] as sequences 
A000079 {n = 3), A090990 (n = 4), A007283 {n = 5), A090992 (n = 6), A000079 
(n = 7), A090994 (n = 8), A020714 (n = 9), A129638 (n = 10), respectively. 
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